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TECHNICAL MEMORANDUM X- 64800 


DOUBLE SCATTERING CORRECTIONS FOR 
THE THEORY OF THE SUN'S AUREOLE 


I. INTRODUCTION 


In our earlier work [1,2] we demonstrated a technique by which we could obtain 
the size distribution (SD), n(r), and the vertical profile, p(y), of the atmospheric aerosols 
from the sky radiance measurements in the region of the sun’s aureole. The theory was 
based on the assumption that only single scattering (SS) is present, which would be the 
case when the sky is clear. However, when the atmosphere in the line of sight of the 
sun is somewhat hazy, the SS approximation may not be valid in general, and higher 
orders of scattering may have to be included to explain the results. 

The precise conditions at which multiple scattering (MS) corrections are needed 
are not yet established. Piaskowska-Fesenkova [3] in an extensive series of experimental 
and theoretical studies has shown that higher order scattering is practically insignificant 
in their aureole results. E. de Bary [4] in her extensive theoretical studies has also 
shown that the influence of multiple scattering is neghgible clossTto the sun. On the other 
hand, according to Van de Hulst [5], SS prevails when the optical depth, r, of the 
medium is less than 0. 1 ; second order scattering or double scattering (DS) should be taken 
into account when 0.1 < r < 0.3; and all orders of MS for r > 0.3. Whether these condi- 
tions should apply to aureole calculations is not clear. To resolve this possible contra- 
diction, we have carried out an explicit calculation of the DS contributions to the radi- 
ance, B 2 , as compared to the SS contribution, Bj . The general approach we take is a 
further development of the work of Dowling and Green [6] and Deepak [7,8] , which are 
rather straightforward but tedious extensions of the first order scattering calculations. 

In the present case we treat the case when the wavelength and atmospheric con- ' 
ii’.’ons arc such that DS is present as a correction. We use the simplifications of the SS 
theory to approximately determine n(r) and p(y) from sky radiance data. Since there 
are certain angular regions in the sky where the DS is negligibly small compared to the 
SS, we make use of this fact in our solution of the problem as explained below. 

The steps for attacking the DS problem are the following: 

1. First, we should get some idea of the angular distance 0° ss from the sun 

m which DS -t SS. This can be done through preliminary computer studies by using 
'fiown models and lest cases of n(r) and p( y). 

2. Next, we interpret the radiance data along the almucantar in this region by 
SS alone (as described in References 1 and 2) and obtain approximations for n(r) and p(y) 



3 . Finally, we substitute these functions back into the expressions for the total 
radiance, Bj + 2 , due to SS and DS, and readjust the parameters in the expressions for 
n(r) and p(y) to best fit the data. 

A brief description of the theory is given in the next section. 


II. THEORETICAL CONSIDERATIONS 


The approach adopted for the DS and the concepts and assumptions used in 
obtaining the expressions for the sky radiance due to the various orders of scattering 
have been explained in detail in References 7 and 8. The basis of our theoretical treat- 
ment is (1 > to treat tire radiance resulting from MS as the sum of contributions due to 

the various orders of scattering, and (2) to consider the Nth order scattering as the result 
of N single scatterings. 

To obtain a clearer physical insight into the problem, we invoke the following 
reasonable simplifying assumptions: 


1 . The atmosphere is composed of air molecules and (various species of) 
aerosols. 

2. The atmosphere is vertically homogeneous; i.e., the altitude-size distribution 
is a separable function of the form 


>?(r,y) = n(r) p(y), [cm' 3 p' 1 ] 


which represents the case when the size distribution, n(r) , of the aerosols does not vary 
with altitude. p(y) is a dimensionless parameter representing the altitude dependence 
ol the total concentration (density) of the particles. 

3. The atmosphere is horizontally homogeneous, which means that in any hori- 
zontal plane at an altitude y , the atmospheric composition remains constant. The com- 
position may change from plane to plane which would be the case of a vertically inhomo- 
geneous atmosphere. Horizontal homogeneity need not imply vertical homogeneity or 


4. Tlie effects of absorption are ignored by choosing to work in the spectral 
regions for which the atmospheric absorption is negligibly small, and only scattering 
predominates. 
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5. The effects of polarization are ignored. This would be a perfectly reasonable 
assumption for SS since we confine our sky radiance measurements to the solar aureole 
region so that it is essentially the case of forward scattering. However, since DS is also 
present, practically all angles of scattering have to be considered. 

6. To obtain the geometrical relations involved, the atmosphere is treated as 
plane parallel, and the effect of the spherical earth, which becomes important whenever 
the scattering angle involved is greater than 75 degrees, is taken into account by the use 
of the concept of generalized Chapman functions [6,9] . 

7. The effects of ground reflection are also ignored in this treatment. 

The effects of absorption and polarization may be introduced into the problem 
by the respective concepts of volume absorption coefficient and the Stokes’ polarization 
parameters. 


111. SKY RADIANCE DUE TO DS 


The geometry of the problem is illustrated in Figure 1. The volume element. dV, , 
which in the case of SS is illuminated by direct sunlight, is now illuminated by light 
scattered by another volume, dV 2 , anywhere else in the surrounding space. The dihedral 
angles between the detector-zenith and zenith dV 2 planes is W] , and that between 
dVj -zenith and zenith-sun planes is co 2 . These are related to the dihedral angle u>, 
between the detector-zenith and the zenith-sun planes by the relation 

CO2 ^ 7T + C0± - U?/ . (2) 


The first scattering takes place at dV 2 and the second at dV, (Fig. 1): 


dV 2 = R 2 2 sin 0 2 d0 2 dec,' dy 2 sec 0 2 


(3a) 


and 


dV, — R, 2 dR,dri - R, 2 sin0, dip 1 dcj| dy : sec0, . (3b 

The scattering angles 0 2 and 0 j , at the dV 2 and dVj respectively, are related to the 
directional angles as given below: 
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cosi// 2 - cosfli cos 02 + sin<?, sin <t> 2 cos (co ( -coj') 


(4a) 


and 



Figure 1 . Geometry illustrating the double scattering 
of sunlight by the atmosphere. 

The altitude dependent volume scattering functions are 


F ' (^2^2) = 


Z F i ^'2 ) 


[km' 1 sr' 1 ] 


and 


(5a) 


F ('P ] ,y i ) Z F j(^i')Pj(y,) , [km’ 1 sr' 1 ] , (: 

where i and j represent the species of scattered that significantly contribute to the 
sky radiance. 
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The optical depths involved are determined for each path just 
analysis. By tedious but straightforward mathematical manipulation 
the total sky radiance due to DS in tha form 


as in the SS 
we may express 


= G(0j,X) / e ’ 2 i T i (y i )D li p'(2) 


(^j.yj)dy] 


( 6 ) 


where 


Djj = - W G = H 0 (X)[e' SiTi(0) W] seC 0, 

F’ (2) (^l,yi> =/ / f F'(^ 2 ,y 2 ) F'f^i'.yj) e 2 i^ T i (y l } ' T i (y 2^ D 2i 


x sec0 2 si n $9 d0 2 du>* dy 2 


(7) 


and 


D 2i " Sj(02 > - Sj(0 1 > 

Sj represents secant functions for angles < 70 degrees, and generalized Chapman 
unctions (defined later) for angles greater than 70 degrees. F< 2 ) may be referred to as 
he equivalent double scattering function” (EDSF). This form, equation (6) is similar 
to the corresponding result for SS, except that the volume scattering function for SS Is 

rep aced b y F ( V^i ,yi ), the scattering function for a unit volume of the atmosphere 
at altitude y, due to DS, which is defined below. To simplify the discussio^e limit 

D and^^Th Scattenng spec,e ® to air molecules and two aerosol components denoted by 
p and q. Thus, in equation (7) we let y 


F 


0A 2 > y 2) - k a( x ) p a^2) p a( y 2) * k p( X)P P ( ^2^p( y 2) 


+ k q( x ) p q(^2>^q( y 2> 


and use a similar expression for F'(^,',y,). 
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The total sky radiance due to SS and DS now may be written 


= B t + B-, = G(0j,X) / e ‘ 2iTi(yi)Dli 

0 

x [F'CiA j ,y j ) + F'( 2 ) i ,y i ) J dy j 


This shows that the presence of DS results in modifying ,y, ) to 


tF'f 0 J ,y 1 ; + F ,f2 )(^|,y ( )l 


Introducing the integrated thickness functions 


F a (D 2i) " w Vq) f e 2 i r i (y 2 )D 2i p a (y 2 )<iy 2 


F p (D a ) ‘ yo> { P p (y2^2 


F <D 9i ) = — 1 — f Q^i T ^ y 2^2i 


V°> o 


q'VT^yi 


where 


w (y) = / P(y) dy , 
y 

and again using tedious but straightforward algebra, we may express B 2 in the form 


B 2 = 2 Bii 
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where i = 1 .2, or 3 denotes a, p, or q, and where each By has the form 



2tt 

f sintf>2 d02 doo', sec^ Jjj 


( 13 ) 


where 

hi = r i (0)r j (0)P i (^ 1 ')P j (^ 2 )F i (D 2i )F(D 21j ) (14) 


and 


D 2 1 j “ D lj ' D 2j 

Here each term yields a DS contribution due to scattering among the various species. 
Let us consider now an almucantar scan (<j > j = 0j), then Djj = 0, D 2 j = hence 

F a (0) = Fp(0) = Fq(0) = 1 and equations (12) and (13) simplify accordingly. 

Notice that when 0, = 0 ls the y-dependence vanishes from the expression for 
the SS radiance, B,, but not from that for the DS radiance. 


IV. SPHERICAL EARTH CQRRECTIONS 


In the computation for B 2 , one encounters the singular secant function since 
integration over elevation angles, 0, from 0 ton is involved, this integrable singularity 
is taken care of by replacing the secants by the generalized Chapman function (GCF) 
(described later in the Appendix and in References 6 and 9) for angles, <f>, close to or 
greater than 7r/2. The GCF is also referred to as the “air mass” in astronomy. It is 
essentially an average of the secant function taken over the size distribution of each of 
the scattering species and is defined by 


/ P(y ( )dy, seoi>' 

S(0,y o ) = 

/ P(yi)dy, 
yo 


(15) 
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where <p is the local zenith angle of the line of sight at volume element P, (Fig i) 

and p(y.) is the density altitude distribution. We have used the detailed approach" 
roughly following that of References 6 and 9 , which considers three regions, <t> < 70 degrees 
< 85 degrees < 90 degrees and a fourth region beyond 90 degrees. 



V. COMPUTATIONAL SCHEME AND RESULTS 


The computations for B 2 involve four integrations 0(0 to n), co(0 to 2rr), 
y,(0 to -) and rj (0 to -). A two-dimensional integration over 0 and w can be per- 
formed very economically and accurately by Conroy’s method [ 10 ]. Integration over 
V, and r, is performed by the use of the 1 2-point Laquerre integration scheme. Com- 

putations were periormed using GCF formulae to take care of the singularity of the sec0' 
tor 0, =90 degrees. 
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We use the oversize distribution function of the form 
N(r) = N 0 [l +(r/a) , ]- i , 


( 16 ) 



where a is a size parameter, v is a dispersion parameter, and N 0 is the magnitude 
parameter. The size distribution is then given by n(r) = - dN/dr. 

The altitude distrbutions are characterized by a form which contains a flat compo- 
nent and a peaked component described by 


h 4 /y-y2>/ h 2 

"i w - TTTFTpq + h ^ + ^ . y2)/h2 j 2 • 07) 

where 0, and 0 2 control the relative importance of the two components, yi is the 
halt-fall-ofl point of the first term, y 2 is the peak position of the second term, and hi 
and h 2 are diffuseness parameters. 

The values of B, and B 2 in the almucantar and sun vertical at various 4/ for 
the August 1 8th event are given in Table 1 . One notices that for the clear day of August 
18th, 1971, the value of B 2 /Bi < 2 percent for 4/ angles up to 20 degrees. The results 
show that in the case of sun vertical scan as one goes toward the horizon (0 L > 0, ) the 
values of B, and B 2 seem to come closer together; i.e., B 2 /B! increases, which is the 
correct trend. The ratio B 2 /B! shows an increase for angles up to 25 degrees as we 
go away from the sun. Our calculations for a less clear day of May, 1970 when t(0) = 

0.6 indicate the DS contribution to the sky radiance in the solar aureole is within 4 
percent of the SS radiance up to angles of 20 degrees. 

Thus, our conclusion tends to support the results of E. de Bary [4] and Piaskowska- 
Fesenkova [3] that multiple scattering is insignificant in aureole studies under reasonably 
clear sky conditions. 
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TABLE I . SS AND DS RADIANCE (THEORETICAL) IN THE 
ALMUCANTAR AND THE SUN VERTICAL FOR EVENT 1 
(AUGUST 18, 1971), GAINESVILLE, FLORIDA 

Data: d ] = 39.6 degrees, X = 550 nm, r a (0) = 0.096, r (0) = 0 182 
r(0) = 0.278 


Almucantar 
Scan (0, = 6 , ) 

Scattering 

Angle 

~~ 

Sky Radiance (Theoretical) 

(W/cm 2 -sr-*im) 

<^1 


B, 

b 2 

b 2 /b, 

(degrees) 

(degrees) 




2 

1.275 

0.667 

l.099(-3) 

0.1651-2) 

6 

3.824 

0.306 

1 .085 (-3) 

0.354(-2) 

10 

6.384 

0.210 

1.062(-3) 

0.505(-2) 

20 

12.71 

0.124 

1.077(-3) 

0.867(-2) 

Sun 





Vertical 





(wi = 0) 





<Pi 





19.6 

20.0 

0.066 

0.659(-3) 

0.99(-2) 

29.6 

10.0 

0.135 

0.826(-3) 

0.61 (-2) 

37.6 

2.0 

0.518 

1.07(-3) 

0.207(-2) | 

41.6 

2.0 

0.543 

1 . 1 7(-3) 

0.21 5(-2) 

45.6 

6.0 

0.243 

1.41 (-3) 

0.579(-2) 

49.6 

10.0 

0.173 

1-57C-3) 

0.910(-2) 

59.6 

20.0 

0.109 

2. 10(-3) 

1.92(-2) 
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APPENDIX 


THE GENERALIZED CHAPMAN FUNCTION 


A brief description and a summary of the expressions for the GCF(s) pertaining to 
the various angular regions are given below. 

If the form of p( y 1 ) in equation (1 5) is exponential, then S(0,y o ) would reduce 
to the usual Chapman function. 

From the geometry (Fig. 2), one obtains 



(A-l ) 


which in terms of the variables 


t = yi - yp 

h i 


and x = 


R c + yp 

hi 


(A-2) 


becomes 


sec 0' = f(t) - 


1 + (t/x) 

[cos 2 <j> + 2(t/x) + (t/x) 2 ] 1 / 2 


(A-3) 


This expression has at 72 singularity near the origin as <S> 1 goes to 90 degrees. 
This is an integrable singularity and can be dealt with by using the forms for the GCF [9] 
that are summarized below. 

For 0 < 70 degrees, the generalized Chapman function 


S (0) = sec <j> 


(A-4) 


For 70 degrees <<j> < 85 degrees, one uses the following expression: 



S(<A, y 0 ) = S e (^, y 0 ) + S r (0, y 0 ) , 


(A-5) 


where, if /3 = cos <p, the essential 


part is given by 


oc 

S e^- yo) = / 

0 


pit) dt 
(j3 2 + 2t)‘ /z 


and the residual part is given by 


S r (0, yo) = f e (t) ^ 



where 




x‘ /z 

(j3 2 + 2t) l/2 


and 


5 t P (t)dt 

O 

t = — 

/ Pit) dt 
0 


For 85 degrees < <j> < 90 degrees, one can use 

oo 

S e (0,y o ) s x J/2 / p(t')dZ , 

0 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


(A- 10) 


where 
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t' = 7a [(Z + 0) 2 - 0 2 ] 

For 90 degrees <0 ^ 120 degrees (Fig. 2b) 

r w Cy b >l (v \ . 

SW.yo) = 2 |_w(yo) J s \ 2 • yt> } ‘ s< ” ‘ * ,y o ) • 

where 


(A- 11) 


(A-12) 


yb = ( R e + yo^ sin ^ - R e (A-l 3) 

and y^ is then greater than zero. 

When y b < 0, the line-of-sight is tangential to or is cut off by the earth’s surface, 
then one uses the form 


S(0, y 0 ) 



S(tt - <f>, yi) - S(jt - 0, y 0 ) 


(A- 14) 


where 

yj = [r 2 2 + 2r 2 (R e + y 0 ) cos 0 + (Rg + y 0 ) 2 ]^-Re . (A-15) 

For 0 > 1 20 degrees, 

S(0 O ) = sec 0 
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